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We consider a model system in which anomalous diffusion is generated by superposition of under- 
lying linear modes with a broad range of relaxation times. In the language of Gaussian polymers, 
our model corresponds to Rouse (Fourier) modes whose friction coefficients scale as wavenumber 
to the power 2 — z. A single (tagged) monomer then executes subdiffusion over a broad range 
of time scales, and its mean square displacement increases as t a with a — 1/z. To demonstrate 
non-trivial aspects of the model, we numerically study the absorption of the tagged particle in one 
dimension near an absorbing boundary or in the interval between two such boundaries. We obtain 
absorption probability densities as a function of time, as well as the position-dependent distribution 
for unabsorbed particles, at several values of a. Each of these properties has features characterized 
by exponents that depend on a. Characteristic distributions found for different values of a have 
similar qualitative features, but are not simply related quantitatively. Comparison of the motion 
of translocation coordinate of a polymer moving through a pore in a membrane with the diffusing 
tagged monomer with identical a also reveals quantitative differences. 

PACS numbers: 05.40.-a 02.50.Ey 87.15.A- 
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I. INTRODUCTION 

In many physical processes one encounters a stochastic 
variable whose mean squared fluctuations increase with 
time t as t a with a / 1. These processes are sometimes 
referred to as anomalous diffusion and specifically 
subdiffusion for a < 1. Such behavior is usually caused 
by the collective dynamics of numerous degrees of free- 
dom, or modes with a broad distribution of character- 
istic times. The exact relations between the underlying 
modes and the observed coordinate are usually unknown, 
and first-principle derivation of the equations governing 
the anomalous diffuser are rare. As a result, a variety of 
such processes are typically grouped into broad classes 
in accordance to their general characteristics @. While 
the exponent a is an important and convenient indicator 
of anomalous dynamics, it contains little information on 
the hidden underlying driving forces. 

A different set of properties of a diffuser can be revealed 
by observing its first passage to a target, or its absorption 
at a trap [3 5]. For instance, it has been established @ 
that the probability density function (PDF) Q(t) to be 
absorbed for one-dimensional (ID) subdiffusion between 
two absorbing walls has a power law tail, if theprocess 
is described by a fractional diffusion equation ||2). This 
slow decay of Q(t) for a subdiffuser following such dy- 
namics leads to a diverging mean absorption time. On 
the other hand, Kantor and Kardar [tJ demonstrated that 
a monomer in a Gaussian polymer that is characterized 
by a — 1/2 has a finite absorption time when it diffuses 
between two absorbing boundaries. The presence of ab- 
sorbing boundaries introduces additional characteristics, 
such as the long time behavior of the survival probability 
S(t) in the presence of a single absorbing wall, or the be- 
havior of the PDF of particle position near the trap. In 



specific cases with a well-defined a (see, e.g., Refs. |8|,l9|) 
some aspects of absorption characteristics have been de- 
termined. However, in the absence of rigorous scaling 
relations one cannot establish whether the exponent a 
determines all other characteristics. Recently, Zoia et 
al. [Io| used some general scaling arguments to propose 
such a relation, furthering the need to probe such prop- 
erties with tunable exponent a. 

While it is generally recognized that underlying (hid- 
den) processes are responsible for anomalous dynamics, 
many theoretical approaches simplify the problem to ef- 
fective equations for the observed variable, hoping to cap- 
ture the multitude of underlying interactions. There is no 
a priori reason for such an approach to succeed, and it is 
therefore useful to consider alternative models where the 
underlying processes are well characterized. In this work 
we consider a model in which subdiffusion is generated 
as a result of superposition of underlying modes with a 
broad range of time scales. For solvability and ease of 
simulation we limit ourselves to linear (but stochastic) 
dynamics for the modes. Our model is closely related to 
the dynamics of a Gaussian polymer, or to fluctuations 
of a Gaussian interface. In the polymer language, which 
we adapt for most of the presentation, the anomalous be- 
havior of a single (tagged) monomer [ll[ is easily under- 
stood in terms of the superposition of underlying Rouse 
modes. The resulting exponent a depends on whether 
the polymer dynamics is diffusive (Rouse) or influenced 
by hydrodynamic interactions (Zimm). The difference 
between the two cases can be cast as due to a wave-length 
dependent friction coefficients. 

In Sec. Ill Bl we take this analogy a step further, and 
show that any value of a can be generated by appropriate 
scaling of the friction coefficients with wavelength. While 
we adapt notions from polymer physics in developing our 
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model, the approach we take is not intended to address 
any particular polymer problem. Rather, we rely on this 
well-defined mathematical model to explore issues perti- 
nent to anomalous diffusion (specifically absorption) , and 
to compare and contrast with other mathematical models 
introduced in this context. Our approach is closely re- 
lated to the model proposed by Krug et al. (Hf where the 
value of a is controlled by modifying the forces between 
particles. We use our generalization to study anomalous 
diffusion in the presence of one and two absorbing bound- 
aries in Sees. IIII Al and IIII B[ In particular, we explore 
the long-time tails of the absorbtion probability Q(t), 
as well as the asymptotic stable shapes of the PDFs of 
the surviving walkers. Qualitatively, various quantities 
have similar features for a variety of a. However, we 
find that a does not enter the results in a trivial way, 
and the stable function for one a cannot be obtained 
from another by simple transformation. Moreover, the 
comparison of results for our walker with exponent a 
coinciding with that of a polymer translocating though 
a membrane pore (Sec. IIV A[) demonstrates quantitative 
differences between the two cases. Some possible exten- 
sions of this work are discussed in Sec. IIVB1 



II. MODEL 

Rouse modes of polymers and anomalous 
dynamics of a monomer 



Polymers I3J provide a relatively simple physical sys- 
tem in which the collective motion of monomers leads 
to behavior spanning a broad range of time scales [Tij . 
Ignoring the interactions between non-covalently-bonded 
monomers, the dynamics of a polymer can be reduced 
to independent Rouse (Fourier) modes [f5j]. For a poly- 
mer consisting of N monomers, each such mode U p 
(p = 0, ••• ,N— 1) has a distinct relaxation time r p . 
The Gaussian (or ideal) polymer is a particularly sim- 
plified model composed of beads (monomers) connected 
by harmonic (Gaussian) springs. Each polymer configu- 
ration is now described by the set of monomer positions 
R n (n = 1 , 2, • • • ,N), and has energy (again neglecting 
further neighbor interactions) 



N 



(1) 



Here the spring constant is k = k B T /b 2 , where k B is the 
Boltzmann constant, T is the temperature, and b is the 
root-mean-square distance between a pair of connected 
monomers. 

From Eq. ([T]) we can construct a simple relaxational 
(Langevin) dynamics for the Gaussian polymer, whereby 

dR 

(- ! ^ = -K(2R n -R n+1 -R n - 1 ) + f n , (2) 
at 

for 1 < n < N. The deterministic force (first term on 
the right hand side) is different for Rx and Rn, since 



the end monomers are attached only to a single neigh- 
bor. Here, £ is the friction coefficient of the monomer, 
while the noise f n (t) has a zero mean and correlations 
of (fn(t)f m (t')) = 2(k B T5 n:7n 5(t ~ t'), to ensure proper 
thermal equilibrium at temperature T. In one dimension 
the positions {R n } are scalars, while the generalization to 
vectorial coordinates in higher dimensions is trivial, as in 
this model the coordinates in different spatial dimensions 
are independent. 

Rouse modes are now obtained by Fourier transforma- 
tion as [lj] 
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N 



Up = jr ^] Rn COS 

71 = 1 



(2n — l)pir 
27V 



(3) 



where p — 0, 1 • • • , N — I is the mode number and U p 
is the mode amplitude. (The choice of cosines automat- 
ically takes care of the modified equations for R\ and 
Rn-) The Rouse coordinates are decoupled and evolve 
according to [3] 



(4) 



where k p = 8 Nk sin 2 and ( p = 2N( for p > 1, 

and £o = N£. The noise W p has again zero mean and 
correlations of 



(W p (t) W p , (tO) = 2k B T( p 5 PtP ,5(t - t') . 



(5) 



Note that the — I internal modes for p ^ behave as 
particles tethered by a harmonic spring, while the center 
of mass (CM), corresponding to p — 0, freely diffuses. 

The linear Eqs. (Q| are readily solved starting from any 
initial condition, and the probabilities for {U p } arc Gaus- 
sians, with a time-dependent mean set by initial condi- 
tions, and a variance 



k R T 



(l-c- 2t/Tp ), for p > I , 



(6) 



dry = 2D CM t . 



The equilibration (relaxation) times of the internal 
modes are t p — Q p /k. p , while the diffusion constant is 



D, 



ksT/N^. There is clearly a hierarchy of relax- 



ation times: the shortest timescale tn-i ~ C/ (4k) is half 
of the time t s = (/(2k) during which a free monomer 
diffuses a mean squared distance between the adjacent 
monomers, b 2 — ksT/n. For the CM, we can define a 
characteristic time tq = b 2 N/D cu — N 2 (/k associated 
with diffusing over the size of the polymer. This is of the 
same order as the longest internal relaxation time t\ ; for 
long polymers To jr\ ~ 7r 2 . 

By inverting Eq. ([3]) one can also follow the position 
of a specific ("tagged") monomer, as 



Rn = U + 2 



JV-1 

E 

P =i 



U p cos 



(2n — l)p7r 
2N 



(7) 
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This equation becomes particularly simple for the central 
monomer c = (N + l)/2 (assuming odd N), as 



(AT-l)/2 

R c = U + 2 ]T (-l) k U 2k . 
k=i 



(8) 



Since each term in the above sum is (independently) 
Gaussian distributed, so is i? c , with variance 



(N-l)/2 

var(i? c ) = cr 2 + 4 cr 2 



2fc ■ 



(9) 



fc=i 



Utilizing Eq. ([S]), we can now distinguish between three 
regimes: 

1. For short times t <C r s , we can expand the expo- 
nential in Eq. © and obtain cr 2 « ftBTt/(N() in- 
dependent of p. The sum in Eq. © then leads to 
var(i? c ) 2D s t, with a monomer diffusion constant 
of D s = fc B T/C. 

2. For very long times t 3> t± all the internal modes 
saturate to a variance that is independent of t. In 
this regime the additional time dependence comes 
from the first term in Eq. Q corresponding to the 
slow diffusion of the center of mass. 

3. The most interesting regime is for intermediate 
times, with t s -C t <C tq, where only the terms 
for which 2t/r2k > 1 will contribute significantly to 
Eq. Focusing on the corresponding modes, we 
obtain 



var(.R c ) « 4 



N/2 



k=k n 



k B T 

«2fc 



dk . 



(10) 



where fc m ; n is determined from the relation T 2 k min — 
2t. Using the usual expression for t p of the 
Gaussian polymer one can immediately see that 
var(i? c ) ~ t 1 / 2 . This interval clearly exhibits 
anomalous diffusion due to the collective behavior 
of the modes; its size can be made arbitrarily long 
by letting N — > oo. 



B. Generalized modes with variable exponent 

In Ref. [7j the tagged monomer was used as a prototype 
of subdiffusion. Of course, as described so far the anoma- 
lous dynamics of the tagged monomer is characterized by 
the exponent a — 1/2. It is possible to modify Eq. (U) in 
different ways so as to produce dynamics for any value 
of a. We shall do so by considering power-law depen- 
dences of the friction coefficients Cp on the mode-index 
p, as motivated by the following two physical models: 

1. Zimm analyzed the motion of a polymer in the 
presence of hydrodynamic flows, which result in in- 
teractions that decay (in three dimensions) as an 



inverse distance between monomers. These inter- 
actions do not change the probability distribution 
in configuration space (which is still governed by 
the equilibrium Boltzmann weight), but do mod- 
ify the relaxation times. Zimm showed [16j that 
the resulting dynamics can be approximately de- 
scribed by Eq. (g]), but with ( p cx p 1 / 2 . We can 
now work through the same steps as in the previ- 
ous section, and find anomalous diffusion for the 
tagged monomer, but with exponent a = 2/3. 

2. Equation (fTJ) can also be regarded as the energy of 
a fluctuating stretched line, with {R n } indicating 
the heights above a baseline |12l[l7j . If the line sep- 
arates two phases of fixed volume, the sum J2 n 
must remain unchanged under the dynamics. Such 
conserved dynamics are mimicked by Eq. ((3]), but 
with ( p oc p~ 2 17]. As long as the noise correla- 
tions satisfy Eq. ([5]) the equilibrium distribution 
remains unchanged, but the dynamics is slowed 
down, such that the fluctuations of a given coor- 
dinate now evolve with exponent a = 1/4. 

Motivated by the above examples, we consider dynam- 
ics according to Eqs. (|4|) and (|5|) with generalized C P given 
by 



C P = 2CiVc(|r) 2 \ forp 



> 1. 



(11) 



Co - CN Z 



The choice of exponent z leads to time-scales 



cap/N? 



c 



4/cshr(p7r/27V) ktH 



(12) 



where the last approximation is valid for p <C N. The 
longest internal mode scales with the number of degrees 
of freedom as N z . This is the conventional notation for 
the dynamic exponent of a fluctuating line, or a free- 
field theory, but somewhat different from that used to 
denote the relaxation of polymers. The Rouse, Zimm, 
and conserved models correspond to z = 2, 3/2 and 4, 
respectively. The dimensionless constant C = 2Ci/Co is 
somewhat arbitrary. It defines the ratio between the time 
characterizing the diffusion of the CM and the relaxation 
time of the slowest internal mode. We chose it in a way 
that for very short times the motion of R c has a diffu- 
sion constant of D s = ksT/C, irrespective of z. For large 
N this leads to C = 1/(2 — 1). We would like to stress 
that Zimm dynamics, as well as other physical systems 
producing anomalous diffusion, are only approximately 
described by Eq. (T4| with length-scale dependent friction 
constants. However, we employ Eqs. (j4|) and (fTTj) as the 
(exact) definition of our mathematical model for anoma- 
lous diffusion with tunable exponents. 

Focusing on the coordinate R c (the tagged monomer) , 
we observe that it executes normal diffusion with diffu- 
sion coefficient D s for times t <C t s = C/2k. For very long 
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times t tq — (£/ ' n)N z it again performs regular diffu- 
sion but with a much smaller center of mass diffusion co- 
efficient of Ds/N 2 ' 1 . At intermediate times r s < t < tq, 
fluctuations of R c are influenced by the p-dependent re- 
laxation times in Eq. (I12p , and evolve anomalously as 
var(i? c ) oc t a . The exponent a can be obtained simply 
by noting [ll| that at times of order To, var(i? c ) should 
be similar to typical equilibrium (squared) size of the 
polymer, which grows as NksT/n. Equating these two 
quantities immediately yields 

f D t\ a 

var(i? c ) - Kb 2 f-^-J , (13) 
where if is a dimensionless pref actor, and 

a=-. (14) 

z 

This result can also be directly obtained from Eqs. (|10p 
and CTJJ). 

In Ref. [12J an alternative strategy is employed for ob- 
taining a tunable exponent, namely by scaling the "spring 
constant" in Eq. 2] clS , while leaving the friction co- 
efficients ( p unchanged. Without a corresponding scaling 
of the noise amplitudes W p , the steady state probability 
is now also modified. The simulations of Ref. [l2[ are 
actually obtained by evolving the Langevin equations in 
real space (as opposed to Fourier mode evolutions per- 
formed in our current work) . This necessitates generaliz- 
ing the interactions in Eq. [2] to further neighbors, and/or 
generating correlated noise in real space. Nevertheless, 
for z = 2 the two approaches should be identical - corre- 
sponding to Gaussian polymers - and direct comparison 
should be possible. 




time 

FIG. 1. (Color online) Mean squared position of the central 
monomer (in units of fa 2 ) as a function of time (in units of r B ) 
for a polymer of length N = 257. The curves correspond (left 
to right) to 2 = 1.25, 1.5, 1.75, 2, and 2.25. Each curve is an 
average over 10,000 realizations. 

Eq. (fT4"|) . Using our model it is possible to get very slow 
dynamics a — > by taking z 3> 1, or to reproduce a = 1 
by setting z = 1. 

We also verified that the probability density of the dis- 
tribution of R c is a Gaussian. The Langevin equation 
describing each U p can be solved analytically, and conse- 
quently the parameters of the Gaussian distribution for 
R c (its mean and variance) are known. Thus, the nu- 
merical results presented in Fig. Q] were anticipated, and 
primarily served to evaluate the accuracy of the numer- 
ical procedure. In the following sections we will use the 
same procedure for results that cannot be found analyt- 
ically. 



C. Numerical implementation 

To verify the dependence of the anomalous exponent 
on z, we simulated the dynamics of a chain of N = 257 
monomers by numerically solving the Langevin equa- 
tion for each of the modes. At the beginning of each 
simulation we equilibrated the polymer by randomizing 
the initial mode amplitudes, and positioned the central 
monomer (c = 129) at the origin. The position R c of the 
central monomer was evolved by numerically integrat- 
ing the Langevin Eqs. (HJ) with the Smart Monte Carlo 
method [llj], followed by transforming the mode ampli- 
tudes Up to the monomer position space, at each time 
step. Figure [T] illustrates the results for several values of 
z. The averages of R 2 were calculated over 10,000 inde- 
pendent simulations. For each value of z the times are 
at least an order of magnitude shorter than the slowest 
relaxation mode of the polymer for that z. For t < r s 
the curves coincide: in that region the particles perform 
normal diffusion with diffusion constant D s which is in- 
dependent of z. For t > t s we can clearly observe a pure 
power law growth with exponent a — 1/z confirming 



III. RESULTS 

The behavior of a particle near absorbing boundaries 
may reveal aspects of the dynamics not apparent in the 
scaling of the mean square displacement. Indeed, con- 
sideration of absorption of a monomer in a simple Gaus- 
sian polymer [7[ have provided novel insights into the dif- 
ferences between different forms of anomalous dynamics 
with a = 1/2. In our numerical studies, we implemented 
anomalous diffusion as described in Sec. Ill C[ but with 
only the tagged monomer interacting with the absorbing 
boundaries. For the case of a single absorbing bound- 
ary, for each z, the starting position x(0) of the tagged 
monomer was at a distance 86 away from the absorbing 
boundary, while for the case of two absorbing boundaries 
the monomer was placed between them at a distance 86 
from either one. The numerical procedure imposes both 
lower and upper limits on the relevant times: 

1. If the tagged monomer is initially located at a dis- 
tance a;(0) from an absorbing boundary, a suffi- 
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ciently long time is required for the probability den- 
sity to be influenced by absorption. Since the typ- 
ical squared distance travelled by anomalous dif- 
fusers is given by Eq. (|13[) . the absorption proba- 
bility becomes significant after a time (disregarding 
the dimensionless prefactor K) 

T = {b 2 /D s )[x{Q)/bf' a = 2r s [x(0)/6] 2 /« . (15) 




10 1 10 2 10 3 10* 

time 



In our simulations with a;(0) = 86, and for z = 1.25 
(2.25) this leads to T w 360r s (2.3 x 10 4 r s ). Al- 
ternatively, one can find directly from Eq. (|9]) , that 
for z = 1.25 (2.25) the quantity {R 2 C ) reaches 64& 2 
at time T « 160r s (2.3 x 10 4 r s ). (Graphically, T 
can be found simply from Fig. [1] as the time corre- 
sponding to 64& 2 .) The fact that T" is not strictly 
proportional to T, is related to the z-dependence 
of K, since z enters Eqs. (jll|l through the constant 
C. 

2. Anomalous diffusion is expected for times signifi- 
cantly shorter than the longest relaxation time tq rs 
t s N z ; for z = 1.25 (2.25) and N = 257, we estimate 
r « 10 3 r s (2.6 x 10 5 r s ). [Since n « r /[(z - 1)tt 2 ] 
(see Eq. (|12p). the corresponding values of the slow- 
est internal modes are 400t s (2.1 x 10 4 r s .)] In the 
presence of a single absorbing boundary our simu- 
lation times for z — 1.25 (2.25) were shorter than 
1.6 x 10 3 r s (1.6 x 10 4 t s ), while for the case of two 
absorbing boundaries they were shorter than 600r s 
(6.5 x 10 3 t s ). These numbers indicate that most of 
our simulations stayed within the anomalous dif- 
fusion regime. To verify this point we performed 
simulations for N = 65, 129, and 257 for z = 2 and 
observed the convergence of their absorption time 
distributions, which indicates that we are in the 
TV- independent regime. In the following we report 
only the results for the central monomer c = 129 of 
the chains with N — 257. 



A. Absorption time distribution 

1. Single absorbing boundary 

The problem of a particle preforming normal diffu- 
sion in the presence of absorbing boundaries has been 
described in detail by Chandrasekhar [19]. It can be cast 
as a simple (linear) diffusion equation for the evolving 
probability density, with vanishing boundary conditions 
at the absorbing points. In the presence of one absorb- 
ing boundary in ID, an elegant solution is found by the 
method of images [j, 0], i.e. by subtracting from the 
Gaussian solution describing the probability density in 
the absence of absorption, a similar Gaussian centered 
at the "mirror image position" with respect to the ab- 
sorbing boundary. At times shorter than T = x 2 (0)/D, 
where D is the diffusion constant of the particle and x(0) 
is its initial distance from the boundary, the particle does 



FIG. 2. (Color online) Logarithmic plot of absorption proba- 
bility distribution as a function of time (in units of r a ) of the 
central monomer of a N = 257 polymer, in the presence of an 
absorbing boundary at a distance 8b from the initial position 
of the monomer. The leftmost curved depicts the result for 
normal diffusion of a single particle. The rest of the curves 
correspond (left-to-right) to z = 1.25, 1.5, 1.75, 2, and 2.25. 
For each z value, 100,000 independent runs were performed. 
The inset depicts the value of 9 obtained from the slopes of 
these graphs as a function of a. The continuous line depicts 
the relation between these exponents proposed in Ref. [2C| 
(see text). 



not feel the absorbing boundary. For t^$>T the survival 
probability (obtained by integrating the solution over the 
allowed interval) scales as S(t) ~ t -1 / 2 , while the absorp- 
tion PDF behaves as Q(t) = -dS/dt ~ t~ 3/2 . We note 
that the mean absorbtion time is infinite, since the par- 
ticle can survive indefinitely by diffusing away from the 
absorbing boundary. 

Figure [2] depicts on a logarithmic scale the PDFs of the 
absorption Q(t) of the tagged monomer initially located 
at distance x(0) = 8b from a single absorbing boundary 
for different values of z. (For comparison, results for nor- 
mal diffusion of a single particle are also shown.) Absorp- 
tion is negligible at short times, but gradually increases 
to a maximum at a time significantly smaller than T de- 
fined by Eq. (Ho} . As in the case of normal diffusion, it 
is generally accepted that the absorption PDF decays as 
a power law Q(t) ~ t^ 1 ^ 6 at long times; 8 = 1/2 for reg- 
ular diffusers while any 8 < 1 leads to a diverging mean 
absorption time. We attempted to extract the exponent 
8 from the slopes of the curves in the logarithmic plot in 
Fig. [5] at large values of t. While 10 4 independent runs 
were performed to obtain each of the curves, only a small 
fraction of diffusers survived to times significantly longer 
than the position of the maximum. Nevertheless, for each 
z we have reasonably accurate results for 1.5 decades be- 
yond the time of maximal Q(t). Only the second half of 
this interval (on a logarithmic scale) is a straight line, 
and it was used to evaluate 8. There are thus significant 
statistical errors in the estimates for 8 as depicted by the 
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error bars in the inset of Fig. [2j 

Studies of continuous time random walks using the 
fractional Fokker-Planck equation [2(| obtain a simple 
relation = a/2 for < a < 2. This relation is depicted 
by the continuous line in the inset. Although there is no 
sound theoretical foundation for applying this relation 
to the collective anomalous diffusion of our model, we 
note that there is some correspondence between the line 
and the measured exponents. In Ref. 12] an alternative 
relation 8 = 1— a/2 is proposed, based on considera- 
tions of fractional Brownian motion. Neither the values, 
nor the trend in this relation are consistent with the re- 
sults in Fig. [2] We are puzzled by this discrepancy as 
Ref. [l2j also provides numerical support for this rela- 
tion, and the method used (although somewhat distinct 
in general) should essentially coincide with ours for the 
case of z = 2 (a = 1/2). We are reluctant to make a 
definite statement regarding this exponent, as in addi- 
tion to the statistical errors (error bars in the inset in 
Fig. [3) there are uncertainties due to possible systematic 
errors: The measurements of 9 are performed for times 
one order of magnitude larger that of the maximum of 
Q(t). These times do not significantly exceed T or T', 
and are possibly even shorter than the latter, especially 
for larger values of z. Getting rid of such possible sys- 
tematic errors requires significantly larger N and more 
statistical samples, but is clearly needed to resolve this 
discrepancy. 



Two absorbing boundaries 



We also considered a tagged monomer confined in the 
interval between two absorbing boundaries separated by 
166. Initially the particle is placed half-way between the 
two boundaries, i.e. at a distance 86 from each. Fig- 
ure [3] illustrates on a semi- logarithmic plot, the absorp- 
tion probability Q(t) for different values of z. The dis- 
tributions have the same general shape as before, with 
absorption probability rising with time to a maximum. 
However, the fall off at long times appears to be expo- 
nential. The straight lines in the semi-logarithmic plots 
clearly rule out the stretched exponential decay charac- 
terizing other forms of anomalous diffusion [2l| . They 
also bear no resemblance to the power law decay (with di- 
verging mean absorption time) expected for subdiffusion 
described by the fractional Fokker-Plank equation [6j. 
The mean absorption time and the typical decay time 
(as determined from exponential decay) are practically 
indistinguishable for each z. They are, however, by more 
than an order of magnitude shorter than the typical time 
T defined by Eq. (15) or V. 
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FIG. 3. (Color online) Semi-logarithmic plot of the absorp- 
tion probability density of a tagged monomer as a function of 
time (in units of r s ) for the central monomer in a N = 257 
polymer, in the interval between two absorbing boundaries at 
a distance 166. The monomer starts at the midpoint between 
the boundaries. Each curve is the result of 100,000 inde- 
pendent simulations. The different curves correspond (left to 
right) to z = 1.25, 1.5, 1.75, 2, and 2.25. 



B. Long-time distribution of particle position 

1. Single absorbing boundary 

Let us now consider the dependence on the coordinate 
x of the PDF p(x,t\x(0)), starting at an initial position 
x(0) from a single absorbing boundary at x = 0. As 
noted before, for a regularly diffusing particle this PDF 
is obtained by the method of images as the difference be- 
tween Gaussians centered at ±a;(0), whose width grows as 
l(t) = y/D s t. Expanding this solution close to the bound- 
ary, we find p(x,t\x(0)) xx(Q)/£(t) 3 , i.e. the PDF van- 
ishes linearly with x in the vicinity of the boundary. It 
is tempting to generalize this result to anomalous diffu- 
sion, by simply replacing the scaling form of the width 
by £(t) = b 1 ~ a (D s t) a / 2 , and again concluding a linear 
behavior with x albeit with a different dependence on t. 
However, the method of images relies on the absence of 
memory in the motion [l^ . which is not correct for our 
non-Markovian processes. Indeed in a previous study we 
observed that for the case of a = 1/2 the vanishing of the 
PDF near an absorbing boundary is faster than linear. 
(A similar problem occurs when a diffuser performs Levy 
flights, since the absorbing boundary is no longer a turn- 
ing point of the trajectory [9, 22].) We shall assume that 
for t 3> T (Eq. (15) ) the behavior near the boundary can 
be described by p(x,t\x(0)) ~ x^. 

In the long time limit we expect £(t) to be the only rel- 
evant length in the problem. However, the initial separa- 
tion, x(0), from the absorbing boundary is another length 
scale, which may become irrelevant only for £(t) ^> x(0). 
To check this assumption we plotted the PDF of unab- 
sorbed particles in terms of the scaled variable p = x/£(t) 
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FIG. 4. (Color online) Probability density function of the cen- 
tral monomer c = 129 in a polymer with N = 257 monomers 
in the presence of one absorbing boundary. The monomer is 
initially located at a distance 86 from the absorbing boundary. 
The horizontal axis is in the scaled variable p = x/£(t). The 
graphs correspond to different times (right to left) t/r s = 2, 
32, 128, 1024, 3500, 6000, and obtained from 100,000 inde- 
pendent runs. 



for N = 257 polymers. Figure [4] depicts a sequence of 
PDFs for z — 2 (a — 1/2) at several times. For short 
times, i.e. when a;(0) S> £(t) the maximum of the PDF 
remains centered close to x(0), and therefore, its center 
appears near p — x(0)/£(t), and moves to smaller values 
of p with increasing £(t). Indeed this process is clearly 
seen in Fig. 2] for the three graphs representing the short 
times, with their maxima moving to the left as £~ 1//4 . 
Graphs corresponding to the two largest times almost 
coincide representing the final stable shape. Note that 
this behavior appears in the same range as the apparent 
power law behavior of Q(t) in Fig. [2] and above T defined 
in Eq. (fT5|), To verify the stability of the scaled PDFs 
it is desirable to study even larger times. Unfortunately, 
the quality of the graphs deteriorates since the number 
of surviving diffusers becomes very small. 

Figure [5] depicts on a logarithmic scale the PDF of the 
particle position in terms of the scaled variable p = x/l(t) 
for several values of z. Since the evaluation of the prob- 
abilities is performed at large times, when only a small 
fraction of the initial 100,000 samples survives, the sta- 
tistical fluctuations are significant. The figure was ob- 
tained by using rather large bins, which poses a problem 
for a function fast approaching zero. In Fig. [3] about five 
leftmost points of each graph containing small numbers 
of events should be disregarded due to statistical uncer- 
tainty and more importantly due to the distortion caused 
by the bin sizes. These effects severely limit the accu- 
racy with which we can determine the exponent <f). As a 
guide to the eye we have added straight lines with slopes 
<f> = l/ct, which seem to provide a fair approximation of 
the slopes in the range 0.1 < p < 1. We note that such a 
form of </> describes the behavior near an absorbing wall 




FIG. 5. (Color online) Probability density function of the 
central particle (c = 129) in chains of N = 257 monomers 
in the presence of one absorbing boundary. The horizontal 
axis is in the scaled variable p = x/£(t) (see text). The dif- 
ferent curves correspond to z = 1.25, 1.5, 1.75, 2, and 2.25 
(bottom to top). The graphs are shifted vertically for clarity 
with increasing z values by a factor of 5. The dashed lines 
have slopes (j> = 1/a. Each curve was obtained from 100,000 
independent runs. 

for Levy flights |8j, as well as for diffusion described by 
a fractional Laplacian Q between two absorbing bound- 
aries. Again, as emphasized before, one should beware 
of possible systematic errors, as the times for which the 
PDF is measured are not particularly long. 

Recently, Zoia et al. argued [l(| that under rather gen- 
eral assumptions there is a relation between the anoma- 
lous diffusion exponent a, the exponent 9 governing the 
tail of the absorption PDF, and the boundary exponent 
0, given by = 20/ a. Relying on = l-a/2 (Ref. flUj). 
they thus obtain <j> = 2/ a — 1 which is larger than the 
estimates from Fig. [5] On the other hand, using our fits 
with 9 w a/2 would yield <j> w 1, which is smaller than 
our data indicates. 



2. Two absorbing boundaries 

In Sec. lIII A~2l we noted that with two absorbing bound- 
aries the absorption probability eventually decays expo- 
nentially. Indeed, the time dependent PDF for a normal 
diffuser between two absorbing boundaries is represented 
by a sum of spatial sinusoidal modes (eigenfunctions of 
the Laplacian operator) multiplied by functions of time 
which are pure exponentials. At large times only the 
lowest harmonic corresponding to the slowest decay, sur- 
vives. Thus for normal diffusion in the interval between 
two absorbing boundaries at Xbi,2 — ±86, the spatial 
probability at long times behaves as ~ cos(7rx/166), again 
vanishing linearly at the endpoints. In Ref. [7j it was 
demonstrated that for a = 1/2, at times significantly 
larger than the mean absorbtion time, the normalized 
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FIG. 6. (Color online) Probability density function of the 
central particle in a chain with N = 257 monomers with ab- 
sorbing boundaries at Xbi,2 = ±86 as a function of position 
measured from the center of the interval (in units of 6) for 
z = 1.25, 1.5, 1.75, 2, and 2.25 (broad to narrow). Each 
graph is the result of 100,000 independent runs. 
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FIG. 7. (Color online) Probability density function of the 
central monomer in a chain with N = 257 particles in the 
presence of absorbing boundaries at Xu = —86 and X\,2 = 
86, for z = 1.25, 1.5, 1.75, 2, and 2.25 (bottom to top). 
Each curve is the result of 100,000 independent runs. The 
curves are shifted vertically for clarity, by a factor of 5 at 
each increasing z. The dashed lines have slopes <fr = 1/a. 



PDF of positions of the surviving tagged monomer has a 
stable shape different from a cosine. 

We performed a detailed study of spatial dependence 
of the PDF of the surviving anomalous walker between 
two absorbing boundaries for several values of z. The 
properly normalized PDF of surviving monomers pro- 
gresses from a Gaussian with variance growing linearly 
in time (for t < r s ), to a Gaussian with variance increas- 
ing as t a at intermediate times, before settling down to 
a stable shape beyond the point of maximum of Q(t) in 
Fig. [3] Figure [5] depicts these stable shapes for several 
values of z. Note the non- linear behavior of the curves 
near the boundaries. Interestingly the boundary expo- 
nent appears to approach <f> = 1 as z — > 1. This is indeed 
the expected behavior for a normal diffuser, although we 
note that our diffusers in the limit z = 1 still reflect the 
collective behavior of many modes. 

To better display the behavior of these stable functions 
near the boundaries, in Fig. [7] we plot them on a logarith- 
mic scale as a function of a distance from one edge. The 
results are distorted not only for the reasons mentioned 
in Sec. IIIIB 11 but also because of the smearing caused 
by the finite time step of the algorithm (causing a typical 
step size of each monomer) . As in the case of a single ab- 
sorbing boundary, the dashed lines with slopes <p = 1/a 
are drawn to guide the eye. Since the curves have sim- 
ilar shapes and approximately follow a power law near 
the boundary, we attempted to collapse various curves 
by raising them to power 2a and normalizing them, but 
this procedure did not result in a good data collapse. 



IV. DISCUSSION 
A. Comparison with translocation 

As we were initially lead to this subject in connection 
with polymer translocation, it is fitting to conclude by 
returning to this issue. Translocation, the passage of a 
polymer through a pore in a membrane, is an impor- 
tant process that has been studied extensively in the last 
decade [2314301 ]. Phages, for example, invade bacteria by 
taking advantage of existing channels in bacterial mem- 
branes to translocate their DNA/RNA inside [3lJ. In 
theoretical models, it is convenient to study a transloca- 
tion coordinate s which denotes the number of monomers 
s on one side of the pore. The dynamics of this coordi- 
nate is anomalous: If we assume that the translocation 
time is of the order of polymer relaxation time To [27| . 
then the variance of s will increase with time as t a with 
a' = 2/{2v + 1) (Rouse dynamics). (For a self-avoiding 
polymer diffusing in 2D, v = 3/4 and a' = 0.8.) We note 
that the actual value of the exponent a' also involves fac- 
tors not explicitly related to the polymer dynamics: (i) 
that self-avoiding effects expand the equilibrium size of 
the polymer in the physical space; and (ii) that the rele- 
vant variable represents a ID coordinate in the internal 
space of monomer numbers. While the expression for the 
actual exponent has been supported numerically in some 
studies [27], H2-H1] , and disputed in others [36-38], the 
anomalous nature of dynamics is not in question. 

To simplify the process, numerical implementations 
frequently begin by inserting half of the polymer into 
the pore (i.e. s(t = 0) ss N/2) and allowing the poly- 
mer to diffuse until either of its ends (s = 1 or s = N) 
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FIG. 8. (Color online) The dots connected by a line depict 
the normalized PDF of the translocation coordinate, at times 
significantly exceeding the mean translocation time, for a two- 
dimensional self-avoiding polymer of length N = 128 that 
did not yet translocate (from Ref. !34j). This is compared 
to the normalized PDF of a central particle (in a chain of 
N = 257) performing anomalous diffusion controlled by ex- 
ponent z = 1.25 (i.e. a = 0.8), moving between two absorbing 
boundaries. For the purpose of comparison the range of the 
translocation coordinate and the range of monomer positions 
have been shifted and rescaled to the segment [0,1]. The in- 
set shows the same quantities on the logarithmic scale; the 
dashed has slope 1.25. 



leaves the pore. This closely resembles the motion of 
anomalous diffuser in Sec. MI B~2l between two absorbing 
walls, and in a previous work [7J for a = 1/2. We are 
now in position to make a more meaningful comparison 
by choosing a value of z that reproduces the observed 
exponent for anomalous dynamics of s(t). Specifically, 
recently Chatelain et al. [3J] preformed high accuracy 
simulations of two-dimensional translocation of a self- 
avoiding polymer. Relevant conclusions from this work 
are: (a) At short times the distribution of s is almost 
an exact Gaussian whose variance increases in time with 
exponent a' sa 0.8. (b) For times larger than the typical 
translocation time the distribution of translocation time 
decays exponentially, as in Q(t) in the current work, (c) 
For times significantly larger than the mean translocation 
time the PDF of the surviving translocation coordinate 
qualitatively resembles those in Fig. [5] of Sec. MI B 21 

In Fig. [5] the results of Ref. [34| are compared with 
simulations of our model with z = 1.25 to reproduce the 
exponent a = 0.8. For better comparison the allowed 
interval in both cases is shifted and rescaled to the range 
between and 1. While the curves are quite similar, 
they do not coincide. The translocation data is repre- 
sented by a narrower bell-shaped curve, and close to the 
boundaries is better described by a power law with ex- 
ponent 4> ~ 1-44 [HI, while the curve obtained in our 
simulations produces a lower exponent of <f> ~ 1.2. The 
differences between the two behaviors is better observed 



on the logarithmic scale in the inset of Fig. [H Thus, 
there are quantitative differences between translocation 
and anomalous diffusion of a monomer with a similar ex- 
ponent a. 



B. Summary 

In this work we concentrated on a group of subdiffu- 
sion processes in which a tunable anomalous exponent 
a is generated through collective behavior of many de- 
grees of freedom. This is achieved by superposition of 
linear modes in which the relaxation times are scaled by 
a power law. In the polymer language this corresponds 
to following a tagged monomer when the friction coeffi- 
cients of the Rouse modes have a power-law dependence 
on wavelength. In the absence of absorbing boundaries 
the model can be solved exactly; starting from a point 
the PDF of the anomalous walker is a Gaussian whose 
width grows in time as t a . We were not able to solve the 
problem in the presence of absorbing boundaries, and re- 
sorted to numerical simulations. With a single absorbing 
point the PDF of absorption decays slowly at long times 
as t^ 1 ^ 8 . The power law decay can be justified by noting 
that the particle can avoid absorption by moving away 
from the trapping point. A qualitative understanding of 
the behavior is not yet attained: Estimates based on the 
fractional Fokker-Planck equation suggest [2(j 9 = a/2, 
while an alternative picture from fractional Brownian 
motion suggests [l2| 9 = 1 — a/2. Reference [l2j provides 
provides numerical support for the latter, while our re- 
sults are more consistent with the former. The possibility 
of systematic errors prevents us from making a definite 
statement on this point, and indicate necessity of further 
work. When the tagged monomer is confined to an in- 
terval bounded by two traps, the survival probability is 
found to decay exponentially at long times, irrespective 
of the subdiffusive exponent. 

An interesting feature of the process is the vanish- 
ing of the PDF on approaching an absorbing boundary 
(whether single or double). The method of images, which 
is only valid for regular (Markovian) diffusion, predicts a 
linear approach to zero, while our simulations indicate a 
singular form characterized by an exponent <fi for anoma- 
lous walks. While we cannot determine this exponent 
precisely due to various sampling problems, our results 
do not appear to support recently proposed exponent re- 
lations The similarities in the shapes of the sta- 
ble PDFs of surviving walkers in an interval for different 
values of a, initially raised the hope that they can be 
collapsed by a simple transformation (e.g. raising them 
to some power). However, the mismatch between the 
curves obtained for different exponents a is sufficient to 
rule out an over-arching super-universality. Furthermore, 
the discrepancies between translocation of a self-avoiding 
polymer and an anomalous diffuser with a similar expo- 
nent, suggest that the exponent a is not sufficient to 
characterize universality. The situation is reminiscent 
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of critical phenomena in which Gaussian (linear) models 
can be devised to reproduce a particular critical expo- 
nent, but which do not capture the full complexity of 
the non-linear theory. The absence of definitive agree- 
ment between numerics and proposed models is on one 
hand disappointing, but on the other hand points to the 
necessity of further work and clarification. 

The linear nature of the underlying model raises the 
hope that exact solutions may be within reach. In the 
meantime the model does provide a means of generating 
anomalous walkers with a tunable exponent that incor- 
porate some realistic features of collective dynamics of 
interacting degrees of freedom. There are certainly puz- 
zles pertaining to the behavior of such anomalous walkers 
close to an absorbing boundary. To answer these ques- 
tions simulations need to probe sufficiently short times 
to remain in the regime of anomalous dynamics, but long 
enough to ensure convergence to stable forms. Our sim- 
ulations had a rather limited range of times satisfying 
the above constraints. While one order of magnitude 



increase in N could open a broad range of validity of 
the above conditions, it would significantly slow down 
the simulations. In addition, working at longer times 
significantly increases the attrition of the samples, and 
requires increasing sample size by several orders of mag- 
nitude. Currently, such ideal conditions are beyond our 
numerical abilities, but some improvement over the cur- 
rent results are certainly possible. 
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